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Asymptotically Safe theories of gravity have recently received much attention. In this work we 
discuss a class of inflationary models derived from quantum-gravity modification of quadratic gravity 
according to the induced scaling around the non-Gaussian fixed point at very high energies. It is 
argued that the presence of a three dimensional ultraviolet critical surface generates operators of 
non-integer power of the type R 2 ~ 0 ^ 2 in the effective Lagrangian, where 9 > 0 is a critical exponent. 
The requirement of a successful inflationary model in agreement with the recent Planck 2015 data 
puts important constraints on the strength of this new type of couplings. 


I. INTRODUCTION 

If the inflationary scenario is correct, then it is possible 
that quantum gravitational phenomena could be detected 
in anisotropy experiments of the cosmic microwave back¬ 
ground radiation and in observations of the large scale 
structure of the Universe [1]. 

Although a consistent quantization of the gravitational 
field is still lacking, it is assumed that quantizing the 
matter fields and the gravitational field on a classical 
background should be sufficient to explain the generation 
of the initial spectrum of perturbations. 

In recent times a promising theory of quantum gravity 
has appeared in the framework of the asymptotic safety 
(AS) proposal. Its central idea, as first discussed by 
Weinberg [2], is to define the continuum limit around 
a non-Gaussian fixed point (NGFP) characterized by an 
ultraviolet (UV) critical manifold of finite dimensionality 
[3-24], 

One of the most striking consequences of this approach 
is the fact that the fundamental theory seems to be rather 
different from the classical gravity based on the Einstein- 
Hilbert action. In the infinite cutoff limit, it is character¬ 
ized by a vanishing Newton’s constant which is therefore 
antiscreened at high energies. Classical gravity can only 
be recovered as a low-energy effective theory at some in¬ 
termediate scale and it does not define a fundamental 
theory. Although we do not know yet the structure of 
the fundamental Lagrangian, all the investigations per¬ 
formed so far, considering general f(R) truncations or 
more complicated tensorial structure like have 

confirmed the finite dimensionality of the UV critical sur¬ 
face (see [25-27] for reviews and also [28] for a recent 
investigation within the f(R) truncation). 

In the simple case of two relevant directions, as in the 
case of quantum Einstein-Hilbert gravity, several investi¬ 
gations have focused on the implications of the running of 
the Newton’s constant in models of the Early Universe. 
In [29-33] it has been shown that the renormalization- 
group induced evolution of the Newton’s constant and 
cosmological constant can provide a consistent cosmic 
history of the Universe from the initial singularity to the 
superaccelerated expansion. 


In particular in [32] it has been argued that the scal¬ 
ing properties of the 2-points correlation function of the 
graviton near the NGFP induce a scale invariant spec¬ 
trum of the primordial perturbations, characterized by a 
spectral index n s which, to a very good approximation, 
must satisfy n s ~ 1. 

Subsequent investigations have tried to produce suc¬ 
cessful models of inflation by considering an extended 
structure of the effective Lagrangian near the NGFP [34], 
but it turns out that the use of an “optimal cutoff” might 
result in a fine tuning problem [35, 36]. In [37] the run¬ 
ning of the gravitational and cosmological constants has 
been described in terms of a Jordan-Brans-Dicke model 
with a vanishing Brans-Dicke parameter, while the via¬ 
bility of the AS scenario in models where the inflaton is 
the Higgs field has been discussed in [38, 39]. 

An effective Lagrangian encoding the leading quan¬ 
tum gravitational effect near the NGFP has been pro¬ 
posed in [40] where a RG-improvement of the linearized 
/3-functions has been performed using the idea that the 
relevant cutoff in this situation is provided by the local 
curvature, i.e. k 2 ~ R (see also [41]). In a recent in¬ 
vestigation a similar approach has been successfully ap¬ 
plied to study the quantum gravity modifications of the 
Starobinsky model [42] (see also [43-48] for other exam¬ 
ples of quantum deformed quadratic gravity inflationary 
models). 

In this work we extend the analysis of [40] by includ¬ 
ing the additional relevant direction produced by the R 2 
operator on the UV critical surface. 

It will be shown that the requirement of a successful 
inflationary model in agreement with the recent Planck 
2015 data release of CMB anisotropy puts significant con¬ 
straints on the renormalized flow generated by the R 2 
term. In particular the dynamics of the inflaton after 
the inflationary phase can be rather different from the 
well known R+ R 2 Starobinsky model. 

The rest of the paper is organized as follows: In Section 
II we introduce the basic formalism and obtain the effec¬ 
tive field theory description in terms of a f(R) theory. 
In Section III we compute the spectral index n s for the 
primordial perturbations and the tensor-to-scalar ratio r 
using a conformal mapping to the Einstein frame. In Sec- 
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tion IV the oscillatory phase after inflation is discussed, 
and Section V is devoted to the conclusions. 


II. BASIC FORMALISM 

Let us consider the quadratic gravity Lagrangian 

L k = -^(R-2X k k 2 )- p k R 2 (1) 

i-Ong k 

where k is a running energy scale and g k , X k and Pk 
are dimensionless running coupling constants whose in¬ 
finite momentum limit is controlled by a NGFP, i.e. 
limfc-xx, (g k ,X k ,P k ) = ( 3 *,A*,/?*) ^ ( 0 , 0 , 0 ) [ 6 , 49]: does 
this running show up at the level of predictions for a 
specific inflationary model? 

In the case at hand the qualitative behavior of the UV 
critical manifold of the quantum theory defined by ( 1 ) 
at the NGFP is rather simple, as its continuum limit 
can be described only by three relevant couplings. In 
particular there exist trajectories that emanate from the 
NGFP and possess a long classical regime where the ef¬ 
fective action is approximated by the standard Einstein- 
Hilbert action [49]. On the contrary the quantitative de¬ 
tails of the renormalized flow around the NGFP are still 
rather uncertain. In fact, not only the precise location 
of the NGFP is regulator dependent (as expected), but 
also the value of the critical exponents depends on the 
truncation strategy employed to solve the flow equation 
[12, 50]. Moreover recent investigations based on uni- 
modular gravity [51], and general arguments [52] suggest 
that the critical exponents are indeed real [53, 54]. 

For these reasons we would like to find an approxima¬ 
tion to the renormalized flow which encodes the general 
qualitative feature of the scaling around the NGFP and 
assumes real critical exponents. We thus approximate 
the running of X k with its tree-level scaling X k ~ cq k~ 2 
where Cq is a dimensionful constant, and decouple the 
running of g k from the running of P k . This latter ap¬ 
proximation is justified by the impressive stability of the 
critical exponents for the Newton’s constant against the 
inclusion of higher order truncations as shown in several 
investigations [ 12 ], but it has also the important advan¬ 
tage to allow for an analitic expression of the flow in our 
model. In fact, under these assumptions the renormal¬ 
ized flow thus reads [ 10 ] 


07rcifc 2 

6ng 2 + 23ci(fc 2 — g 2 ) 

(2) 

®3 

/ 7.2 \ 2 

Pk = P* + bo \ ) 

( 3 ) 


the critical exponents for the R 2 coupling. 

It is important to stress that, as long as C\ < 67 t /23 
the running described by Eq. (2) smoothly interpolates 
between the Gaussian fixed point (GFP) and the NGFP, 
g* = 9 k{k —> 00) = 67 t /23 and therefore it captures the 
qualititive features of the flow described in [49] . The con¬ 
stants 60 , Co and c\ depend on the physical situation at 
hand and must be fixed by confronting the model with 
observations: in principle these are the only free param¬ 
eters of our theory corresponding to the three relevant 
directions of the UV critical surface of the action (1). 
In other words, by changing ci, Co and bo it is possi¬ 
ble to explore various RG trajectories all ending at the 
NGFP. The relevant question is if it possible to actually 
constraints these numbers, in particular the value of bo. 


As argued in [40, 42, 55] by substituing ( 2 ), (3) and 
X k = co/k 2 in ( 1 ), a renormalization group improved 
effective Lagrangian can be obtained by the scale identi¬ 
fication k 2 —> £R, where £ is positive number. It must be 
stressed that the general structure of the resulting RG- 
improved effective Lagrangian agrees very well with the 
high-curvature solution of the fixed point equation for 
a generic asymptotically safe f{R) theory, as it emerges 
from the analysis of [20] and [40]. 

At last, we obtain the following RG-improved action 



n s 3 

R + aR 2 2 


R 2 

6 m 2 



( 4 ) 


where g is chosen so that k 2 = 8ttGn = 487r 2 c\ / (67r g 2 — 
23 (g 2 + 2£co)ci) = 1/M p i 2 . Moreover A = g 2 c o(67r — 
23 ci)/(67r/i 2 — 23 (g 2 + 2£co)ci) and M p 2 /to 2 = 
12(23£/(967t 2 — /?»)); in particular if £ > 967r 2 /3*/23 and 
Co < g 2 (6n — 23ci)/46£ci, then 1/m 2 and A are positive 
definite. 


On the other hand 

a = -2// 3 & 0 /M pl 2 (5) 

which only depends on 0 3 and bo, but not on the fixed 
point values. Concerning 0 3 the numerical evidence ac¬ 
cumulated so far has shown that its value is rather stable 
against the introduction of higher order truncation in the 
flow equation [50], as it should be expected for a critical 
exponent. On the contrary bo is by construction a non- 
universal quantity whose value cannot be determined by 
the RG group. It labels a specific trajectory emanating 
from the fixed point and its actual value should be de¬ 
termined by matching with a low energy observable. 

Let us now introduce an auxiliary field ip defined via 


where g is an infrared renormalization point, c 3 = g k {k = 
g) (see [10] for details). According to [ 6 ] /3» = p k (k —> 
00) ~ 0 . 002 , while b 0 is a free parameter obtained by the 
linearization of the RG flow around the NGFP, and 0 3 is 


«’<*> Sl + a ( 2 — + (6) 

In principle this relation can always be inverted at least 
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locally, provided ip iX ^ 0 so that (4) is equivalent to 



ifR -{ip- 1) x(p) 


—a 




x(<p) 

3 m 2 


(7) 


In practice, due to its non-linearity, the task of inverting 
Eq. (6) can be very difficult and one must often resort to 
numerical methods. In our case, according to the analy¬ 
sis of [12], as 63 is rather close to unity we can safely set 
03 = 1 for any practical calculation (as it can be numer¬ 
ically checked). In this case we explicitly obtain the two 
branches 



(' TJc?m 1 + 8 m 2 <p — 8 m 


2 


± 3^3 yj27a A m s + 16 a 2 m 6 {ip — 1)) 


( 8 ) 


with the reality condition x > 1 — 27m 2 a 2 /16. 

Using these solutions we can obtain a canonically cou¬ 
pled scalar field by introducing a conformally related 

metric g pv by means of g^ = pg^v with p = eV 2 / 3 ^. 
In the Einstein frame 1 [57, 58] the action thus reads 


r e - ^ 9 e V d^d v (l> - V±{<t>) 

(9) 

where 

m 2 e -2 \/i K 0 f / /y , \ 2 

V±W= 256 k 2 { 192 ( eV ^ X ) - 3a + 128A 


S = J d^X\J—gE 


-3a 2 (a 2 + 16 e^ K,p - Hi) =F 6a 3 \J a 2 + 16 e^^ - 16 


-V32a 


a* + 8eV - 8 ) ± aV a 2 + 16eV f *<t> - 16 


( 10 ) 


and we have measured a and A in units of the scalaron 
mass m by means of the rescaling a a/3V3 m and A —► 
Am 2 , so that both a and A are dimensionless numbers. 


III. INFLATION DYNAMICS AND 
PRIMORDIAL SPECTRUM 

For A 3> 1 the potential behaves as V±{<p) ~ 
1 e^ 2 V^ K ^m 2 A/k 2 while for a 1 we have V±(<j>) ~ 


1 Note that, as discussed by [56] the field redefinition employed 
in going from the Jordan to the Einstein frame involves new 
additional contributions in the path-integral arising from the Ja- 
cobians. In investigations including the presence of matter field 
in the starting Lagrangian, these new terms cannot be neglected. 



Figure 1. Inflation potential for various values a and A = 1.4. 


±a 4 e 2 V^ k ^to 2 /k 2 . In both cases this implies as(f) ~ 
t 3 / 8 which does not provide inflation. 

In order to study the inflation scenario in slow-roll ap¬ 
proximation, we need to know the shape of the poten¬ 
tials V±{(p). First we notice that, for all values (a, A), 
the potential has the plateau V±{<f>) = for </>—>• oo. 
To verify slow-roll conditions inflation must start from 
(f> > M p i and then proceed from the right to left. The 
behavior of the potential for (f> <C M p i, and thus the in¬ 
flation scenario, strongly depends on the values (a, A). 
In particular V±{<p) can either develop a minimum, or 
V±{4>) — > —oo for <j) —> —oo. In this work we study the 
class of potentials such that the inflation ends after a 
finite number N of e-folds, with a phase of parametric 
oscillations of the field <j>. In order to have a well defined 
reheating phase it is clear that the potential must have 
a minimum; furthermore it can be proved that, in our 
case, a well defined exit from inflation occurs only for 
potentials with a mininum cf> m j n such that V(<^ m i n ) < 0. 
These conditions are verified only for V + (<p) if a £ [1,3] 
and A £ [0,1.5] (the potential V-{(f>) can have a mini¬ 
mum for some special values of (a, A), but it is always 
V{<t>m i„) >0). In these cases the potential is depicted in 
fig. (1) for various values of a. 

In other words, although for a and A very close to zero 
(10) is only a small modification of the classical Starobin- 
sky R + R 2 model, for a £ [1,3] and A £ [0,1.5] the 
potential V+{(j>) develops a non-trivial minimum at neg¬ 
ative values of the potential which makes our model sig¬ 
nificantly different from the original R + R 2 Starobinsky 
inflation. In particular, as we shall see, after exit from 
inflation, the dynamics of the preheating phase is charac¬ 
terized by a lower limit to \(j)\ as it can be seen in fig. (3) in 
the </>-</> plane (see [59] for a general study of inflationary 
potentials with a negative value of the minimum). 

Let us introduce the slow-roll parameters e and 77, 


e(0) 


J_ (V^)\ 

2« 2 l V{</>) J 


V{4>) = 7A 


1 (v"{4>) 


V v{<t>) 


(ii) 


and define the number of e-folds before the end of infla- 
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Figure 2. We compare the theoretical predictions in the r- 
n 3 plane for different values of a for the Planck collaboration 
2015 data release for the TT correlation assuming ACDM + r 
[60]. Triangles are for N = 55 and squares for IV = 60 e-folds. 
Solid and dashed lines are the 1 <t and 2 a confidence levels, 
respectively. 


are displayed in the following table: 


Cases 

N 

= 50 

N 

= 55 

N 

= 60 

A 

a 

n s 

r 

n s 

r 

n s 

r 


1.0 

0.965 

0.0069 

0.968 

0.0058 

0.971 

0.0050 

0 

1.8 

0.966 

0.0074 

0.969 

0.0063 

0.972 

0.0055 


2.6 

0.967 

0.0076 

0.969 

0.0065 

0.972 

0.0056 


1.0 

0.965 

0.0070 

0.968 

0.0059 

0.971 

0.0051 

1 

1.8 

0.966 

0.0074 

0.969 

0.0063 

0.972 

0.0055 


2.6 

0.967 

0.0076 

0.969 

0.0065 

0.972 

0.0056 


Confronting n s and r with the T5 Planck data, our model 
agrees very well, as displayed in fig. (2). 

According to Eq. (13) the normalization of the scalar 
power spectrum at the pivot scale fc* = 0.05 Mpc -1 pro¬ 
vides us with m ~ (1.5 T 7) • 10 14 GeV, depending on the 
value of a. 

It should be stressed that, in our model, if A < 0 there 
would be no exit from inflation with the standard reheat¬ 
ing phase. As argued in [61] the presence of matter field 
could change the sign of A depending on the number of 
Dirac, scalar and vector fields. We hope to address this 
issue in a future investigation. 


tion 

C^end 

N= H(t) dt. (12) 

J tN 

The amplitude of the primordial scalar power spectrum 
reads 


IV. OSCILLATORY PHASE AFTER INFLATION 

After the end of inflation, the inflaton field (f> begins to 
oscillate around the minimum 0 m j„ of V+(<J)). To study 
this phase, we can do the following approximation 


2 

K 247r 2 e((f>) 


(13) 


and, in this approximation, the spectral index n s and the 
tensor-to-scalar ratio r are given by 


V+tt) ~ V(</>) = | [(</> - 0 mi n) 2 - b] (15) 

Where </> min (a, A), a(a,A) = P" (<?W) and b(a, A) = 
—2 V+(^ min )/y"(0 m in) depend on the values of a and A, 
and are given by the following expressions 


n s = 1 — 6 + 2 r](cf>i) 

r T = 16 e((f>i) 


•y/| (3a 3 (a 2 — 4) — 32aA + 4 (a 2 — 6) |a| 3 ) 
0 min (a, A) = 6a (a 2 _ 8 ) ( a 2 + 2 ) - 64aA + 8 (a 2 - 9) |a| 3 


Where <f>i = is the value of the inflaton field cj>(t) 

at the beginning of the inflation. 

In our case it is not possible to analytically evaluate 
these expressions, and we must resort to a numerical es¬ 
timation of the initial value of the field which determines 
the final number of e-folds at the end of inflation. As 
usual inflation is supposed to stop when the slow-roll con¬ 
dition is violated, e(4>) = 1. This condition determines 
the final value </>/ = <fi{t e nd) of the inflaton field, while 
the initial value of the field (f>(t) is obtained by inverting 
Eq. (12), once a number of e-folds is fixed. Our results 


. .. 48 + 18a 2 — 3a 4 + 32 A — 4a 3 |a| + 36a |a| 

“(a.A) = - — -^-— 

8a (15a 4 - 3a 6 - 96A + 8a 2 (15 + 4A)) |a| 
_ | (48 + 18a 2 - 3a 4 + 32A - 4a (a 2 - 9) |a|) 2 
-25a 8 + 132a 6 - 384a 2 A 
| (48 + 18a 2 - 3a 4 + 32A - 4a (a 2 - 9) |a|) 2 
48a 4 (21 + 4A) - 1024A(3 + A) 

| (48 + 18a 2 - 3a 4 + 32A - 4a (a 2 - 9) |a|) 2 
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d,<f>( t) Mb) 



r , and A = 10. This case arise when V+((p m in) > 0. 

1 igure 3. Phase space evolution m the (p-cp space tor a = 1 

and A = 1. Note the presence of a limit cycle. 


The time evolution of the field is given by the Fried¬ 
mann equation 

4>{t)+3H(t) <j>(t) + V'{</>(£)) =0 (16) 


Where H(t) is the Hubble constant 


3 H(t) = 3 




1/2 


The scale factor a(t), in the case ab = 0 (P(</> m in) = 0, 
Starobinsky model), describes the usual matter dom¬ 
inated era, while the solutions with V(</> m i n ) 7 ^ 0 are 
compatible with a matter dominated era only at the 
beginning of the oscillatory phase. On the other hand 
a consistent treatment of the following reheating phase 
must include the contribution of the matter fields, an 
investigation which is beyond the scope of this work. 


\j\ [<P(t) 2 + a (</>(t) - ^ min ) 2 - ab] 1/_ (17) 


V. CONCLUSIONS 


Putting x{t) = y'a (</>(t) — </> m in) and y(t) = 4>{t), equation 
(16) reads 


y = - [§ (y 2 + x2 - ab )] 2 V - V a x 

x = yfay 


(18) 


The long time behavior of this dynamical system is 
mainly determined by the sign of ab = —2 U(<^ min ). If 
ab < 0 (i.e. V (</; m j rl ) > 0) then the point (0, 0), that is the 
minimum point (<j> m i n , V^min)), is an attractive node. If 
ab > 0 (i.e. P(0 m in) < 0) a limit cycle ( y 2 + x 2 = ab) 
appears. Note that, in this sense, ab = 0 is an Hopf bi¬ 
furcation point. A numerical study of the original Fried¬ 
man equation (16) has confirmed our analytical findings. 
In particular we have obtained the phase diagrams rel¬ 
ative to the cases U + (0 m m) < 0 (our case, fig. (3)) and 
V + ((j) m in) > 0 (fig. (4)). The dynamical system analysis 
can be useful to determine the scale factor a(t) after the 
inflation era, that is related to the Hubble constant by 
H[i ) = In particular, using the polar coordinates 

and then performing the method of averaging, we obtain 


sin \ab\ tj 
a(t) = <j t 2 / 3 

sinh ( \/| |afe| t 


2/3 


2/3 


ab > 0 
ab = 0 
ab <0 


(19) 


In this work we have extended the idea of [40] by 
including in the renormalized flow the presence of 
the additional relevant direction associated to the 
R 2 operator. We have approximated the flow of this 
operator around the NGFP with its linear expression, 
an approximation which should capture the essential 
qualitative features of the flow and allow an analytical 
treatment of the resulting non-linear f(R) Lagrangian. 

The most important point of our investigation is 
that our inflation model should significatly differ from 
the well known Starobinsky model because it predicts 
a tensor-to-scalar ratio which is significantly higher, 
and an inflationary dynamics which is characterized 
by a limit-cycle behavior at the inflation exit. More 
important, our predictions are in agreement with the 
latest Planck data which put important constraints on 
our model. 

An important limitation of this study is the simple ten- 
sorial structure of the effective Lagrangian which assumes 
a functional dependence of the f(R) type. Quadratic op¬ 
erators like R^R^ are also associated to relevant di¬ 
rections around the NGFP and in principle their pres¬ 
ence could dramatically change the inflation dynamics 
and the generation of the primordial spectrum of fluc¬ 
tuations. On the other hand, our approach can easily 
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be extended in order to include the contribution of these 
additional operators, whose critical exponents have been 
calculated in [13]. We hope to address this issue in a 
future work. 
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